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ABSTRACT 

The distance of the Fermi-detected blazar 7-ray emission site from the supermassive black hole is a 
matter of active debate. Here we present a method for testing if the GeV emission of powerful blazars 
is produced within the sub-pc scale broad line region (BLR) or farther out in the pc-scale molecular 
torus (MT) environment. If the GeV emission takes place within the BLR, the inverse Compton 
(IC) scattering of the BLR ultraviolet (UV) seed photons that produces the 7-rays takes place at the 
onset of the Klein-Nishina regime. This causes the electron cooling time to become practically energy 
independent and the variation of the 7-ray emission to be almost achromatic. If on the other hand 
the 7-ray emission is produced farther out in the pc-scale MT, the IC scattering of the infrared (IR) 
MT seed photons that produces the 7-rays takes place in the Thomson regime, resulting to energy- 
dependent electron cooling times, manifested as faster cooling times for higher Fermi energies. We 
demonstrate these characteristics and discuss the applicability and limitations of our method. 
Subject headings: galaxies: active — quasars: general — radiation mechanisms: non-thermal — 
gamma rays: galaxies 



1. INTRODUCTION 

Blazars dominat e the extragalactic 7-ray sky 
(|Abdo et al.l l2011aD . exhibiting flares in which the 
Fermi 7-ray flux is see n to increase up to several times 
the pre-flare level (e.g. lAbdo et al.l[2010b( ). Because the 
emitting regions are unresolved, the location of the GeV 
emission remains unknown and this constitutes an open 
issue in understanding energy dissipation in blazars. 

In the leptonic model, the GeV emission of blazars is 
produced by IC scattering of photons off of the same 
relativistic electrons in the jet t hat giv e rise to syn- 
chrotron radiation (for a review see iBottc hcr 2007). The 
seed photons for IC scattering are synchrotron photons 
originating within the relativistic j et (synchrotron - self 
Compton, SSC scattering) (e.g. IMaraschi et al.|[l992D or 
photons originating external to the jet (external Comp- 
ton, EC scattering), such as UV accretion disk photons 
(iDermer et al.' IT992l). UV photons originating from the 
BLR (Sikora e t al.l 1199^ . and IR photons originating 
from molecular torus (MT) (e.g. lBiazeiowski et al.ll2000( ). 

Because the 7-ray emission of blazars can vary sig- 
nificantly on short times, down to the instrument lim- 
its, ranging from a fe w minutes (in the TeV range) 
dlh aronian eFaH [2009h to a few hours (in the GeV 
range) (p^'oschini et al.ll20lTl ). light-time travel arguments 
constrain th e size of the emitting region and it has been 
argued (e.g. iTavecchio et al.|[2010D that this requires the 
emission to take place at the sub-pc scale. Explain- 
ing the Fermi GeV spect ral breaks seen in powerful 
blazars (jAbdo et al.ll20103 ) as a result of pair absorp- 
tion on the He Lyman recombination continuum and 
lines, produced in the inner part of the BLR, places the 
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GeV emission well within the BLR ([Poutanen fc SternI 
l2010f ). Very short distances < 10^^ cm), however, where 
accretion disk photons dominate are not favored, be- 
cause of the high pair production 7-ray opacity (e.g. 
iGhiseUini fc Tavecch io 2009). 

Gamma-ray flares have been associated with radio and 
optical flares, placing the GeV e mission outside the BLR: 
i Lahteenmaki fc Valtaolal (|2003[ ). using the temporal sep- 
aration of the radio and 7-ray fluxes, argued that the 
7-ray emission site is at a distance of ~ 5 pc from the 
central engine. Connections have been also drawn be- 
tween 7-ray outbursts, optical flare polarization changes, 
and ejections of superl uminal radio knots on the pc scale: 
iMarscher et al.l (|2010f ) associated short 7-ray variations 
with the passage of disturbances from the radio core, 
parsecs away from t he black hole. In the case of OJ 287, 
lAgudo et al.l ()2011D found 7-ray flaring to be simulta- 
neous with mm outbursts associated with a VLBA jet 
component 14 pc downstream from the radio core. 

Here we propose a method to distinguish if the GeV 
emission of powerful blazars takes place within the sub- 
pc BLR or further out within the pc-scale MT. In fj2] 
we discuss our understanding of the BLR, MT, and SSC 
photon fields and where each one dominates. We show 
that cooling on the UV BLR photons results in achro- 
matic cooling times in the Fermi band, while cooling 
on the MT IR photons results to faster cooling times at 
higher 7-ray energies. In f|3]we present our method, dis- 
cussing the dilution of an energy-dependent cooling time 
from light travel time effects and discuss the applicability 
of our method. Finally, in ^l^e present our conclusions. 

2. COOLING OF GEV EMITTING ELECTRONS 
2.1. Characterizing the BLR and the MT 
Reverberation mapping of the BLR points to a typi- 
cal size of Rblr. ~ 10^^ ^^145 cm (e.g. iKaspi et al.ll2007l : 
iBentz et al.l 120091 ). where Ld.45 is the accretion disk lu- 
minosity in units of 10''^ erg s~^. Because of this scaling, 
and assuming that Lblr = isLRLd, where (,blr ~ 0.1 
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is the BLR covering factor, the energy density inside the 
BLR {Ublr ^ 10~^ erg cm~'^) i s not expected to vary 
widely from object to object ('Ghiseh ini fc Tavecchiol 
[2009). As was shown by iTavecc hio (200^ when trans- 
formed to the jet frame, the BLR spectrum, dominated 
by H Lya photons, is broadened and boosted in lu- 
minosity due to relativistic effects and can be well- 
approximated by a blackbody spectrum that peaks at 
vblr = l-hTvLya {^BLR « 3 X 10"^ F), whcrc F is the 
bulk Lorentz factor and e is the photon energy in mc^ 
units. 

I n the MT the dust temperature ranges from T ^ 300 
K (ILandt et all [20lol) I' to T ~ 1200 K (IClearv et al.l 
I2007t iMalmrose et al.ll2011f) . with hotter emitting ma- 
terial expected to b e closer to the central engine (e.g. 
[Nenkova ct al. I2008D . Because of the larger size of 
the MT, it has been probed by reverberation mapping 
only for Seyfert galaxies (e.g. ISuganuma et al.l 120061) 
which have lower luminosities and therefore smaller MT 
sizes. These st udies, along with near -IR interferomet- 
ric studies fe.g. iKishimoto et al.ll20Hl ) also support an 
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scaling, suggesting that U is similar for 



sources of different luminosities. Adopting the scaling 
of IGhisellini fc Tavecchiol (l2009l ). Rmt = 2.5 x IQ^^l]/^^ 
cm and assuming a covering factor ^ 0.1 — 0.5, 

the MT photon energy density is Umt ^ 10^'' erg 
cm~^. We also adopt a temperature Tmt ~ 1000, corre- 
sponding to a blackbody peak at vmt ~ 6.0 x 10^^ Hz 
(eMT « 5 X 10-7). 



2.2. Which seed photons dominate where? 

For an isotropic photon field, the co-moving (jet frame) 
energy density U' is U' « (4/3)F^[/, while for pho- 
tons enterin g the emitting region from behind, U' ~ 
(3/4)F-2[/ (jPermer fc Schlickeiseri flOOl . If the emis- 
sion site is located at i? < Rblr, the BLR photon 
field can be considered isotropic in the galaxy frame 
and, using the BLR energy density discussed in ^^i^ 
UUt.r 1-3 X F^o erg cm'^. Similarly, the MT pho- 



BLR 

ton field is isotropic inside the BLR and its co-moving 
seed photon energy density is C/ji/j- ^ 1.3 x 10~^ F^q 
erg cm^'^. Clearly, at Rblr scales, U'^j^j^ dominates 
over U'j^jrp by a factor of ~ 100. If the emission site 
is located at i? ~ Rmt, then the BLR photons en- 
ter the emitting region practically from behind, so that 
U'blr ~ 7.5 X IQ-'^T^o erg cm'^. The MT photons re- 
tain the same co-moving energy density as before and, 
in this case, the MT dominates the co-moving photon 
energy density. 

These external photon field co-moving energy densities 
need to be compared to the synchrotron one. If Rbiob is 
the size of the emitting region, then the co-moving syn- 
chrotron photon energy density is U'g ~ Lg/ {AttcRIi^i^S'^) , 
where d is the usual Doppler factor and Lg the syn- 
chrotron luminosity. An upper limit to Rbiob is set by 
the variability timescale t^: Rbiob ^ ctyd/{l + z), where 
z is the source redshift. We then obtain a lower limit 
U', « L,{1 + 0)V(47rc3t2<56) j// _ g_3 ^ lO'^L^ 45(1 + 

z)H 

v\d^w '^'^ where ty^id is the variability time 



in days. Setting (5 = F the condition U'^ < U'^^f. requires 

•3L,(1 + Z)2y oe(LsAd^ + zYY 



F > 



8.6 



u' 
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This translates to F > 4.8 x Ll^^^^^il + zy/%^/f for 

1 /s 

the emission site at Rblr scales and F > 8.6 x Lj^f^{l + 

z)^/*t~^^'* for the emission site at Rmt scales. In general, 
powerful qua sars exhibit superluminal motions requ iring 
F - 10 - 40 (jJorstad et al.ll2005t iKharb et al.l[2010l ) and 
therefore, the dominant source of seed photons for these 
sources should be the BLR or the MT, provided the 7- 
ray emission site is at a distance not signifi cantly larger 
than t he size of the MT. A recent finding bv lMever et al.l 
(|2012f) that the ratio of 7-ray to synchrotron luminos- 
ity increases with in creasing radio core, as expected fo r 
EC emission ([Derme r 1995; Georganopoulos et al.ll200il) . 
supports further the possibility that powerful blazars are 
EC emitters. 

2.3. Cooling in the BLR vs Cooling in the MT 

The critical difference between the BLR and the MT 
is the energy of the seed photons: the BLR produces 
UV photons while the MT produces IR photons. This 
difference by a factor of ~ 100 in typical photon energy 
is critical in that it affects the energy regime in which the 
GeV-emitting electron IC-cooling takes place, and thus 
the energy dependence of the electron cooling time. 

For this consideration to be relevant, IC cooling must 
dominate over synchrotron cooling. This condition is sat- 
isfied, since in powerful blazars the IC luminosity clearly 
dominates over the synchrotron luminosity, reaching in 
high states an IC luminosity higher t hat the synchrotron 
one by a factor o f up to ~ 100 (e.g. lAbdo et"aril2010al : 
IMever et alllMl) . 

For electrons cooling in the Thomson regime (769 << 
1, where both the electron Lorentz factor 7 and the seed 
photon energy eo are measured in the same frame) the 
cooling rate 7 cx 7^. For electrons with 760 >> 1 cool- 
ing takes plac e in the Klein -Nishina (KN) regime with 
7 oc ln7 (Blu menthal fc Gou ld 1970). For the broad in- 
termediate regime, 10~^ < jeo < 10^, a parametric ap - 
proximation has been suggested (jModerski et al.l [20051) . 
In what follows we calculate numerically the electron en- 
ergy loss rate 7, which for monoenergetic seed photons 
of number density no and dimcnsionless energy cq is 

j^no r^f{x){e-eo)de, (2) 
./n 4en 



where, following iJonesI ()1968D . 

f{x) = 2a:: log X X 1 - 2x'^ 



(4eo7a:) 
1 -I- 4eo7a; 



(1-x) (3) 



with X = e/[4j'^eo (1 - e/7)]. 

The effects of the transition between Thomson and KN 
regimes on the electron energy distribution (EED) and 
the resultant spectrum of the synchrotron an d IC emis- 
sion have ^ been studied before (e. ^. Blumen thal 119711: 



Zdziarski' '1989^; 'Dermer fc Atovan"2002; Sokolov et alj 
2004; .Moderski et al...2005; .Kusunosc fc Takahara. 20051 



Localizing energy dissipation in blazars 



3 



Thomson Cooling 

Klein-Nishina Cooling 

Synchrotron Cooling 




10' 10' 10* 



Y 



Fig. 1. — Electron cooling time in the galaxy frame as a function 
of 7. Top panel: Blazar emission site located in the BLR. Bottom 
panel: Blazar located in the MT. Broken lines represent the various 
cooling mechanisms and the solid black line is the total cooling 
time. Plots were calculated for seed photon energies eg BLR = 
3 X 10~^ and cq^mt = 5 x 10"'^, seed photon energy densities 
Ublr = 10"'^ erg cm~^, Umt = 10"* erg cm""^ and magnetic 
field energy Ub = ^OOUec- 

GeorganoDOulos et alJ 120061 : iManolakou et alJ 120071 : 
Sikora et al.ll2009li ! 

In short, the transition from Thomson to KN cool- 
ing in the EED can be seen for the steady state 
71(7) DC / (5(7)^7/7: for a power law injection (5(7) oc 
7~P, the cooled EED hardens from 71.(7) oc 7^(^+1) 
in the Thomson regime to n{-/) oc j~^p~^'> in the KN 
regime. The transition takes place gradually at 7 ~ 
eQ^. A second transition back to n(7) cx 7~(p+i) is 
expected at higher 7, as synchrotron losses (assumed 
to be less significant than Thomson EC losses) be- 
come larger than KN losses. The imprint of these con- 
siderations on the observed synchrotron and IC spec- 
tra has been considered for the case of the X-ray 
emission of large scale jets (Dermer & Atovan" '2002*) 
and in the case of blazars j Kusunosc & Takahara 2005; 
IGeorganopoulos et al.ll2006t ISikora et al.il2009l ). 

Spectral signatures, however, are not unique and can 
be diluted or altered by extraneous causes (e.g. the 
accelerated electron distribution has intrinsic features 
that deviate from a power law, the synchrotron emission 
is contaminated by the big blue bump, the IC emission 
is modified by pair production absorption from syn- 



chrotron or external photons). For these reasons their 
use in understanding where the GeV emission takes 
place is problematic. 

3. LOCATING THE GEV EMISSION SITE 

There is however another observable aspect of the 
Thomson - KN transition that is free of the issues faced 
by relying on spectral signatures: the energy dependence 
of the electron cooUng time can be used to evaluate the 
regime at which the electrons that produce the GeV radi- 
ation cool, and through this evaluate where the emission 
takes place. 

In the Thomson regime, the electron cooling time 
tc — 7/7 scales as 7"^, while in the KN regime the 
electron cooling time scales as as 7/ In 7. These two 
regimes connect smoothly forming a wide valley around 
760 ^ 1 in which the cooling time is practically energy- 
independent. The energy of the seed photons originating 
in the BLR is greater than that of photons from the MT 
by a factor of ^ 100. As a result, (see Fig. [T|) the lo- 
cation of the energy-independent valley of the BLR case 
is manifested at electron energies lower by a factor of 
^ 100 than that of the MT case. In turn, the energy 
dependence of tc causes a related energy dependence in 
the observed falling time tf of the flare. The falling time 
tf of Fermi light curves can be used, therefore, to deter- 
mine whether a flare occurs within the BLR or within 
the MT. 

3.1. Flare decay in the BLR vs flare decay in the MT 

We utilize an one-zone code to demonstrate the ef- 
fect of a flare occurring within the BLR versus a flare 
occurring within the MT. The code assumes the injec- 
tion of a power law electron distribution and follows 
its evolution taking into account the KN cross section, 
SSC losses, and assuming a black body external photon 
field. T he implicit numerical scheme is sim flar to that 
used b y IChiaberee fc Ghisellinil dTool) and IGraff et al.l 
(|2008[) . After the code reaches a steady-state, an in- 
crease in the injection rate is introduced that lasts for 
a short time before it returns to its initial value. The 
code was initialized to simulate the flare occurring within 
the BLR and the MT respective, described by the val- 
ues discussed in ^ We assumed a Compton dominance 
L I c, max /L synch, max « 50, sourcc size i? = 3 X lO^'^cm, 
bulk Lorentz factor F = 10. Deep in the Thomson 
regime (MT case) the Compton dominance corresponds 
to UecJUb, where Uec is the external energy density 
and Ub is the magnetic energy density. However, as 
cooling transitions from the Thomson to the KN regime 
(BLR case), the ratio Uec/Ub over-estimates the ob- 
served Compton dominance by a factor of ^ 2. 

In the BLR case (Fig. [2J blue curve) IC cooling oc- 
curs at the onset on the KN regime, resulting in a near- 
flattening of t/ within the Fermi band, where a change 
in energy by a factor of 100 (from 200 MeV to 20 GeV) 
results in a change in cooling time of less than an hour, 
with all the Fermi energies having practically indistin- 
guishable decay times (between ~ 2 — 3 hours). In the 
MT case (Fig. [21 red curve) cooling occurs mostly in 
the Thomson regime and t/ is heavily energy-dependent; 
the same change in energy by a factor of 100 results in 
a change in cooling time by ~ 10 hours. This is also 
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Fig. 2. — 1/e falling time in the galaxy frame versus observing 
energy (GeV) for MT seed photons (red), BLR seed photons (blue) 

apparent in the light curves shown in Fig |3l A flare 
occurring within the BLR (Fig [3J top panel) has compa- 
rable values oitf within the Fermi energy range. A flare 
occurring within the MT exhibits distinct falling times 
that decrease with observing energy ([Sj lower panel). 

For suiEciently bright Hares with short decay times 
(comparable to the cooling times of few to several hours 
anticipated), Fermi light curves can be generated at 
multiple energies, and the energy dependence (or lack 
thereof) of tf can be used to reveal the location of the 
GeV emitting site. 
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Fig. 3.— Light curves at 200 MeV, 2 GeV, and 20 GeV for BLR 
seed photons (upper panel) and MT seed photons (lower panel), 
normalized for a steady state L=0 and maximum L=l (in arbitrary 
units). Times are in the galaxy frame. The dashed black line is 
the 1/e decay time. Times are in the galaxy frame. 



3.2. Light Travel Time Effects 

Because the GeV emitting region is not a point source, 
the observed 1/e decay time tobs of the light curve is 
a convolution oi tf and the blob light-crossing time tic- 
Any other variations (e.g. a gradual decrease of the elec- 
tron density) will have a similar effect. For BLR achro- 
matic cooling, we expect the light curves to remain achro- 
matic after considering light crossing time effects. To 
see the effects of tic in the case of MT cooling, where, 
tf cx e~^/^, we consider a source with exponentially de- 
caying emission coefficient, j(t) = jo exp~*/*-f , and calcu- 
late the expected light curves for two energies differing 
by 10 and for a range of tic- As can be seen in Fig I3.2[ 
although tobs increases with increasing tic, the difference 
of tobs between cle and she is practically preserved for 
all tic- Moreover, because for given e^E and enE the ra- 
tio ottf is also given, observations of the 1/e decay times 
for eLE and she can be used to find both tic and tf for 
each energy. 

3.3. The feasibility of the diagnostic 

Very br ight flares, such as the November 2010 flare of 
3C 454.3 ()Abdo et al.ll2011bf) have enough photon statis- 
tics to produce light curves in two Fermi energy bands. 
The practical application of our diagnostic requires that 
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Fig. 4. — Exponentially decaying light curves of a flare in the 
MT. High and low energy luminosities {euE I^LE = 10) a^r^ plotted 
(normalized to 1) for increasing ti^. Solid lines represent euE and 
dashed lines e^f;. Colors represent ti^ at differing factors of t^, 
where tf^HE is normalized to 1. The dashed black line is the l/e 
decay time. 
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the decay times in different energy bands have statis- 
ticaUy meaningful differences if the flare occurs within 
the MT. To evaluate if this is the case for a flare as 
bright as that of 3C 454.3, we adopted the maximum 
observ ed errors in the 0.1 — 1 GeV range and > 1 GeV 
range (lAbdo et al.ll20lfbh and produced simulated data 
assuming energy-dependent cooling in the form of expo- 
nentially decaying light curves. We fitted an exponen- 
tial decay function to the simulated data for each energy 
band, with the characteristic decay time as a free pa- 
rameter. The decay times we recovered were statistically 
distinguishable, suggesting that our diagnostic is appli- 
cable, at least for the brightest flares observed. 

3.4. Constraints from upper limits of the decay time 
difference 

For flares from which no statistically significant differ- 
ence in the decay time can be established, constraints 
can still be imposed on the location of the Fermi de- 
tected emission by considering Atmax , the maximum de- 
cay time difference allowed by the data, between a high 
energy enE and a low energy eLE- The requirement that 
Atmax > tfi^LE) — tfi^HE) cau bc castcd as upper limit 
for the location R of the Fermi emission 

2 1 

3%/3 7rmec2(l + z)i/2ei/^(e-]^/2 _ J 

(4) 

where L mt is the MT luminosity. In terms of quantities 
used in i i2.1l and for eLE, £le corresponding to 100 MeV 



and 1 GeV, i? < 2.3 x lO^^ rio{Atmax,hLMT,45/{'^ + 
zY^^y/^cm. The Lorentz factor can be estimated from 
VLBI studies (e.g. lJorstad et al.ll2005[ ) and the MT lumi- 
nosity can be taken to be a fraction ^~ 0.1 — 0.5 of the 
accretion disk l uminosity or in some cases can be directl y 
measur ed (e.g. lLandt et al.ll2010l : fMalmrose et al.ll2011| ). 

4. CONCLUSIONS 

We presented a diagnostic test for determining if the 
GeV emission in powerfull blazars originates from within 
or outside the BLR. Our method utilizes the fact that if 
electrons cool via IC scattering on BLR photons, cooling 
occurs at the onset of the KN regime and the resulting 
cooling times (and thus light curves) should be achro- 
matic. If cooling occurs via IC scattering of less energetic 
MT photons, the electron cooling times, and thus light- 
curves, should exhibit significant energy-dependence and 
decrease as photon energy increases. The energy depen- 
dence of the cooling time difference perseveres in the 
presence of energy independent timescales such as the 
light crossing time. The method is applicable for the 
brightest Fermi 7-ray flares, and can provide upper lim- 
its on the flare location even when the decay times in 
different energies are statistically indistinguishable. 
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